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ABSTRACT 


A numerical and experimental investigation of the interaction of a 
pair of vortices with a free surface has been undertaken. The analysis 
is based on the vortex-sheet representation of the free surface and the 
use of the appropriate boundary conditions. The experiments were 
performed in a large basin and the vortices were generated through 
the use of a special nozzle. The rise of the resulting Kelvin oval, the 
trajectories of the vortices, and the instantaneous shape of the free 
surface were recorded on a video tape and then carefully analyzed 
through the use of a Motion Analysis system. The results have shown 
that the rise of the vortices not only gives rise to two scars, with a 
pronounced hump in the middle, but also, and more importantly, to a 
three-dimensional instability heretofore unknown. The measured and 
calculated vortex trajectories and the free-surface shapes at the corre- 
sponding times and Froude numbers are found to be in reasonable 
agreement. The new instability will form the basis of future 


investigations. 
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I. INTRODUCTION 


This is part of continuing basic research toward the understand- 
ing of the fundamental mechanisms and physical processes underlying 
two- and three-dimensional vortex/free-surface interactions in homo- 
geneous, stratified, and sheared media, taking into account ambient 
turbulence, viscous effects, and various large-scale instabilities (sinus- 
oidal instability and vortex breakdown) for ship- and submarine- 
related hydrodynamics in a real ocean environment. 

The name “remote sensing” was first used in the early 1960s, 
when it was recognized that the term “aerial photography” was no 
longer suitable to describe new kinds of images formed using energy 
outside the visible spectrum. The developments in 1970s (largely by 
NASA) gave rise to the development of a serious scientific interest in 
remote sensing. This interest led to the routine availability of multi- 
spectral data of large portions of the earth’s surface. As far as the 
oceans are concerned, the interpretation and use of the multi-spectral 
data required a clear understanding of the physics of various phenom- 
ena leading to surface signatures resulting from the body-wake free- 
surface interaction in two- and three-dimensional, sheared, and/or 
stratified media, with or without ambient turbulence. 

As far as the ship wakes are concerned, the most significant wake 


signature has been thought, for nearly a century, to be the Kelvin wave 


pattern, with its 38.9° wake angle (see Figure 1). Observations made 
with remote sensing devices (e.g., synthetic aperture radar) have 
changed this perception dramatically and have shown clearly that a 
narrow strip of disturbed wake (scars and striations) exists for several 
hours behind the ship (see Figures 2a and 2b) with little or no sign of 
the Kelvin wake. These features have been termed turbulent wakes, 
although it is not clear that turbulence alone is responsible for their 
appearance. While the understanding of the physics of this phe- 
nomenon, through experimentation and numerical simulation, has 
attracted a great deal of attention, a Satisfactory model describing the 
essential features of the radar images has not yet been presented. It 
appears that many competing and interacting processes are responsi- 
ble for the observed scar and striation pattern. It has been suggested 
that vortices produced by the ship’s hull may be responsible for the 
suppression of surface waves near the ship track and an enhancement 
of the waves near the edges of the smoothed area. Not only the free- 
surface signatures generated by the vortical flow but also the natural 
waves radiate energy away from the wedge-shaped wake region and 
thereby affect the flow itself. It appears that the dark, narrow images 
are visible due to a decrease in radar return (negative spectral pertur- 
bation) aft of a ship and not to an increase in radar backscatter along 
the edges of the wake. 

The three-dimensional temporal development of the surface sig- 


natures resulting from the interaction of an ascending inclined vortex 











500 m 


Figure 2a. SAR Image of the Narrow Turbulent Wake of a Ship 








100 m 


Figure 2b. Close-Up View of the SAR Image of the Narrow 
Turbulent Wake of a Ship 


pair with the free surface was first described by Sarpkaya [Ref. 1] and 
Sarpkaya and Henderson [Ref. 2]. When a trailing vortex pair generated 
by a lifting surface approaches the free surface, with or without mutual 
induction instability and/or vortex breakdown, the vortices and/or the 
crude vortex rings give rise to surface disturbances (scars and stria- 
tions, see Figures 3a and 3b). The striations are essentially three- 
dimensional free-surface disturbances, normal to the direction of the 
motion of the lifting surface, and come into existence when the vortex 
pair is at a distance equal to about one initial vortex separation from 
the free surface. The scars are small, free-surface depressions on 
either side of a nearly symmetrical hump, directly above the plane of 
symmetry of the vortex pair. They come into existence towards the 
end of the pure striation phase and when the vortices are at a distance 
equal to about sixty percent of the initial vortex separation (b,)) from 
the free surface. The scars are created and driven by the vortices and 
both move outward. Sarpkaya [Ref. 1] and Sarpkaya and Henderson 
[Ref. 2] found that the vortices diffuse rapidly and their circulation 
decreases slowly at first (prior to the occurrence of scars) and rapidly 
thereafter due to the counter-sign vorticity generated in the overlap- 
ping regions of the vortex pair in the recirculation cell and at the free 
surface. 

Previous studies of a vortex (or vortex pair) in ground- or free- 
surface effect have been generally concerned with the boundary layer 


development, generation of oppositely-signed vorticity, the flow 
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resulting from a single vortex held fixed in a uniform flow, and the so- 
called “rebounding” phenomenon, according to which the vortices 
move at first toward the boundary, as on an inviscid trajectory, and 
then away from it [Refs. 3-6]. The boundary surface may be either a 
rigid boundary at which the no-slip condition is satisfied or a clean 
free surface corresponding to zero shear stress. In either case, the 
boundary surface gives rise to oppositely-signed vorticity and the net 
circulation of the vortices decreases with time [Ref. 5]. 

Sarpkaya and Henderson's [Refs. 1, 2] theoretical analysis of the 
scars created by the trailing vortices was based on the classical solu- 
tion of Lamb [Ref. 7], assuming the vortices to be two-dimensional and 
the free surface to be a rigid plane. For small Froude numbers F 
(F = Vos/gbo, where Vo is the initial mutual induction velocity of the 
vortices and bg is the initial vortex separation), each scar front was 
shown to coincide with the stagnation point on the Kelvin oval, formed 
by one of the pair of the trailing vortices and its image. For Froude 
numbers larger than about 0.15, not only the deformation of the free 
surface but also the nonlinear interaction between the said deforma- 
tion and the motion of the vortices are significant. Thus, the free sur- 
face may no longer be assumed rigid. Even though it was fully realized 
at the outset that the problem to be ultimately solved is the under- 
standing of the effects of three-dimensionality, shear, stratification, 
and the parameters characterizing the body-shape and motion on the 
evolution of the surface signatures, the relative ease of the study of the 


two-dimensional flow attracted the immediate attention of 


experimentalists and numerical analysists alike in the past few years 
(Sarpkaya, et al. [Ref. 8]; Willmarth, et al. [Ref. 9]; Dahm, et al. [Ref. 
10]; Marcus [Ref. 11]; and Telste [Ref. 12]). 

Sarpkaya, et al.’s two-dimensional experiments were conducted in 
a 12-foot-long, three-foot-wide, and four-foot-deep (12' x 3' x 4') water 
basin. The nearly two-dimensional vortex pair was generated through 
the use of two counter-rotating plates. Neutrally buoyant fluorescent 
dye was used to visualize the flow. Sarpkaya, et al. [Ref. 8] also carried 
out a numerical analysis of the scars through the use of vortex dynam- 
ics. Their results have shown that the two-dimensional evolution of the 
free surface can be calculated (up to certain normalized times) 
through the use of line vortices or dipoles. The analysis becomes more 
robust as the Froude number is increased. At larger times, the results 
become unrealistic because the laminar and turbulent diffusion of vor- 
ticity are ignored. No wave train was observed on either side of the 
scars for any of the Froude numbers encountered in the analysis or 
experiments. 

Willmarth, et al. [Ref. 9] and Dahm, et al. [Ref. 10], at the Univer- 
sity of Michigan, adopted Sarpkaya, et al.’s [Ref. 8] idea of counter- 
rotating plates to generate vortices and carried out similar 
experiments on vortex-pair free-surface interaction. They too have 
used vortex methods to calculate the path of the vortices and the 
deformation of the free surface. However, no detailed comparisons 
were presented between the observed and calculated scar 


characteristics. 
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Marcus [Ref. 11], working with Berger at UC Berkeley, used the 
finite-difference method to investigate the inviscid, two-dimensional 
interaction between a pair of counter-rotating line vortices and a free 
surface. He has not conducted any experiments and encountered con- 
siderable numerical-stability problems at lower Froude numbers. At 
higher Froude numbers, his predictions were somewhat larger than 
those measured and calculated by SarpkKaya, et al. [Ref. 8]. Telste 
[Ref. 12] used vortex dynamics to investigate the same phenomenon 
(inviscid two-dimensional interaction between a two-dimensional vor- 
tex pair and the free surface). At medium and lower Froude numbers, 
his calculations required extensive numerical filtering to avoid or to 
delay numerical instability. The calculations are generally terminated 
when the free surface forms a sharp corner or when difficulties are 
encountered in obtaining rapid convergence. There does not Seem to 
be any direct relation between the breakdown of the numerical model, 
breakdown of the waves, and the attainment of maximum scar condi- 
tions. Telste offered no comparisons with experiments. 

It is clear from the foregoing that the investigation first initiated 
by Sarpkaya [Ref. 1] at the Naval Postgraduate School quickly mush- 
roomed into a full-scale scientific effort on vortex/free-surface 
interactions, partly because of its intrinsic interest and partly because 
of its far-reaching practical applications in ship- and submarine- 
related hydrodynamics in a real ocean environment. It is also clear 


that much remains to be accomplished, both theoretically and 
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experimentally, toward the understanding of the evolution of surface 
signatures. 
The present investigation has a two-fold purpose: 

1. To conduct further parametric studies of the variables affecting 
the performance of the numerical model through the use of the 
vortex-sheet representation of the free surface; and 

2. To investigate in greater detail the motion of a vortex pair 


approaching and intruding through the free surface and to delin- 
eate the characteristics of the scars and striations. 


PZ 


II. NUMERICAL SIMULATION 


A. DESCRIPTION OF THE MODEL 

The flow model consists of two counter-rotating vortices (CRVs) 
of constant circulation, moving towards the free surface, and a vortex 
sheet, represented by a number of equi-spaced line vortices at the free 
surface. The model is a generalized boundary integral/vortex method 
developed by Baker, Meiron and Orszag [Ref. 13]. 


The use of Bernoulli’s equation on either side of the fluid interface 


leads to 

oa s(u2, + u2,) + gn, + _ = F(t) (1a) 
and 

B24 Hut, +u2,) + gmat B2 =F, (1b) 


in which 96 is the velocity potential, u is velocity, g is gravitational 
acceleration, t is time, n is Surface elevation above mean water level, p 
is density, p is pressure, s is distance along surface, the suffixes t and 
n denote components tangential and normal to the surface, and the 
suffixes 1 and 2 denote the lower and upper surface. Furthermore, on 
the surface one has 

P= Po =P 


Uni = Un2 = Un 
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Tis Naa 
% = Ut1 -uL 
ut = (uty +ut2)/2 (2) 
= (uj, + ug)/2 
V = (v) + v9a)/2 
I =y, ds = ydx 


The use of the foregoing plus the free surface condition given by 


De 
Dr Ont oy) +F n=0 (3) 


lead [Ref. 14] to the momentum equation for the calculation of the 
evolution of vorticity 


O 
8 Dery s raf eeu (Mv vWjm 
C O O 
+ 2 csy2 | 5%] 


(4) 


where the Atwood number At = (92 - 91)/(p2 + 91). For the air/water 
interface, At = -l. 

The velocity of a vortex may be decomposed into two parts: One 
stemming from the use of line vortices (i.e., the two CRVs and the 
surface vortices, equi-spaced at mid-points of the segments) and the 
other from the corrections to the first due to the distributed nature of 


vorticity along the sheet. Thus one has: 


Uk/Vo = (U1 + U2)/Vo (Sa) 
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Vi/Vo = (V1 + V2)/Vo (Sb) 


where the U; and V; components of velocity are obtained from the 
complex function w(z). Normalizing the strengths of all vortices by the 
strength of the counter-rotating vortex pair, all velocities by 


Vo (= To/2nbo), and all distances by bo, one has 


w= “In(z -Z.)- ~In(z +Z)- LT Intz =A) LT ,Inz 1 a) (8) 


k=] k=] 


in which Z) and -Zp, are the positions of the CRVs and z, and -Z, 
denote the positions of the free surface vortices. The velocities U; and 


V, are then given by 


(U) -i Vi)/Vo = dw/dz (7) 


and the velocity components U2 and V2 are given by [Ref. 14] 


ax{ dy an 121) 8 
v= S(t , = 
Ax ( on 
oar Tie 2 sai 


k 


From a mathematical point of view, one would like to have the 
free surface extend from -<o to te and the CRVs originate at (0.5, ~o) 
and at (-0.5, -c). This is impossible from a numerical point of view. It 
is necessary to obtain guidance from experiments to limit the com- 


puter time. Observations have shown that the free surface rises in a 
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relatively small region directly above the CRVs and no waves are gen- 
erated beyond the scars (see Figure 4). In fact, the linearized analysis 
of Kochin, et al. [Ref. 15] of the motion of a line vortex under a free 
surface has shown that the resulting steady wave train far from the 
vortices would be negligibly small. The sample calculations led to the 
conclusion that the free surface can be restricted to a region extend- 
ing from x = -l0 to x = +10. Furthermore, the CRVs are assumed to 
originate at a depth of y = -5, in view of the fact that the free surface 
does not begin to deform until the CRVs arrive at about y = -1. 
Theoretically, one can calculate the potential function 6 from the 
complex velocity potential w, the velocity of the vortices from dw/dz = 
u-iv (excluding the vortex whose velocity is being calculated), and the 
elevation of the free surface from equations 3 or 4 for a given Froude 
number, and trace the evolution of the free surface as a function of 
time. This relatively simple-sounding procedure is anything but sim- 
ple, primarily because the numerical short-wavelength instabilities 
eventually lead to large-scale instabilities. The critical time for the 
onset of these instabilities decreases with decreasing Froude number. 
It is because of this reason that the preliminary calculations were ini- 
tially restricted to large Froude numbers. To be sure, the instability is 
not limited to boundary-integral methods. It is also an inherent char- 
acteristic of the space discretization methods. In fact, the choice of a 
particular method depends, more or less, on the user’s past experi- 


ence, the availability of a computer code, the objectives of the 
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Figure 4. Photograph of a Typical Free-Surface Deformation 
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investigation, and the particular problem. A lucid discussion of the 
advantages and disadvantages of various techniques is presented by 
Yeung [Ref. 16]. Previous investigators dealing with marching methods 
for nonlinear water wave problems have used various filtering [Ref. 17] 
or regridding techniques to eliminate or delay the instability. In gen- 
eral, it is not easy to assess the accuracy of the results. Sensitivity 
analysis, parametric studies, and comparison with experiments may 
help to assess the error introduced by numerical smoothing and fil- 
tering techniques [Ref. 18]. One must also bear in mind that the 
numerical model might not necessarily correspond to a physically 
realizable flow (turbulent diffusion, axial instabilities, wall effects, etc.). 
In the present investigation, it was decided to use Rosenhead’s 
[Ref. 19] smoothing technique, according to which any complex veloc- 
ity is, in general, desingularized by introducing a cut-off distance 64, 


given by (see also Sarpkaya [Ref. 18)]): 


u.-iv, ==) -——*>- ———__- (9) 


The use of other smoothing techniques (e.g., the Longuet-Higgins 
and Cokelet technique [Ref. 17]) was rejected because it tends to over- 
smoothen the free surface. Furthermore, as argued by Roberts 
[Ref. 20], the five-point scheme of Longuet-Higgins and Cokelet may 
lead to hopelessly large errors in some applications. According to 
equation 6, the solution should approach the inviscid flow solution as 


56 —+— 0. Numerous calculations have shown that a 6/Ax value between 1.1 
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and 1.5 yields nearly identical results. The final calculations were 
performed using 6/Ax = 1.25. 

The specific details of the numerical calculations are as follows: (i) 
calculate the strength of the vortex sheet and the strength of the sur- 
face vortices by iteration until dy/dt and/or dU/dt and dV/dt converge to 
an acceptable limit, (ii) find the velocity of all vortices (including those 
of CRVs) and convect them for a time interval At using a second order 
modified Eulerian convection scheme (At = 0.01 in the calculations), 
(iii) calculate the area between the free surface and the x axis, and (iv) 
repeat the above steps. The procedure described above is relatively 
simple and does not require excessive computer time (about 30 min- 
utes on IBM 3033). However, the free surface eventually does become 
unstable. The instability shows up as a kink in the free surface caused 
by the propensity of the neighboring vortices to orbit about each other. 
Some of the instabilities can be alleviated without smoothing. Initial 
calculations have shown that the vortices near the y axis begin to slide 
sideways as the free surface (or the vortex sheet) stretches. Conse- 
quently, the thinly populated regions of the sheet do not yield a 
smooth surface and the flow begins to leak between the vortices. Also, 
the regions where the free surface is depressed (where the scars 
develop) become overpopulated with vortices, leading to the growth of 
short-wavelength Helmholtz instability. This problem can easily be 
alleviated either by packing the vortices more densely near the y axis 
at the start of the calculations, so that the entire surface becomes 


more or less uniformly represented at later times, or by rediscretizing 


Jee 


the vortex sheet at Suitable time intervals. This latter technique has 
been used in the present calculations (for additional details see 


Sarpkaya and Shoaff [Ref. 21)). 


B. RESULTS OF THE NUMERICAL CALCULATIONS 

Numerous runs were made at various Froude numbers from 0.35 
to 0.65. In each case, the calculations were carried out until the insta- 
bilities were strong enough to cause a numerical floating point over- 
flow error at some point during the execution of the program. A 
numerical filter was not used to eliminate or to delay the instabilities. 
Nondimensional time was used as well aS nondimensional lengths and 
speeds. For the sample results shown here, the counter-rotating vor- 
tices were placed at dg/bDo = 3. The start of their motion, measured by 
T = (Vot/byo - do/bo), was T = -3. The free surface was not allowed to 
start to deform until the vortex pair was within 2.5 units of the free 
surface. Horizontal placement of the vortex pair was at -0.5 and 0.5 
units so that the distance between the two would be bo. 

Figures 5 through 8 show the evolution of the free surface at rep- 
resentative times for a Froude number of 0.53. The corresponding 
velocity-vector plots are shown in Figures 9 through 12. The path of 
the main vortex follows the path of a vortex pair approaching a rigid 
boundary until about T = -1.6. Then it starts to diverge and, instead of 
turning to the right as in the case of the rigid wall, it starts to move to 
the left very slowly. This is because the strengths of the discrete vor- 


tices placed near the center line are inducing a horizontal velocity 
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Figure 5. Evolution of Free Surface at T = 0.95 
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Figure 6. Evolution of Free Surface at T = 0.85 
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Figure 7. Evolution of Free Surface at T = 0.75 
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Figure 8. Evolution of Free Surface at T = 0.65 
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Figure 9. Velocity Field at T 
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Figure 10. Velocity Field at T 
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Figure 11. Velocity Field at T = 0.75 


27 


Y/BO 


Vas uf, 
A AG 
AAA LR 
444 \ ENS LN 
444 WMO NN NV VN VV VV V VV VV 
ZV SS 
oy VV KV OVO VOY ey vey ayy 
Van yyy YY VVVVYVVV VV YY VV VV 
i Sa a a co 2 2 2 2 0 
OE inh yy hey uy 
Nee eae YY eee 
NN we eee LAE LEE PP ee 
BN yee ec ec ec 424 4Lebebbeerererrrvrrrrrereryp 
no i a a a 2 2 2 2 2 2 2 A 
AR RPeeec ec ec ce eZ e 4 LA LEE bebe bebe eae ae ba 
ARR Yew e cece ee 4 eA LEE KE KEELE EEE EL 
ARKPRYrxr ce ce Ke eS LL LL EEE EEE 
ARRKXeYT STK KK KK SS SL 4 4LELLELEL 
ARR YY eves ce cee Fe eZ LLL 4 4E4LELEL 
ARKK¥eTT TC KKK KK SEES EZ LEESE 
ARNRRPKYTTT KK KK KEK FH 4 EEE EES 
ARR RS SYS ES © CS Hg Be Bee 
ARAKKRPNEHTTTTST SKS SK KL LS SS SS 
ARAKRHBSPRTSTTS TSK SF EEK HEF LEAFS 
A RN KAN S SS} } SS SS See 
ARRNKRHFeYENNTWT TS SS SK FF FFF FF 
ARN RN WR ® & YS SS SS Sm ee Se eK 
ANAANNRANKRBRPRHEKTHETTTN TS TS SHEE EES 


DO 
CA 
pes 
4 


Figure 12. Velocity Field at T = 0.65 
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component to the left. This eventually would lead to the vortex moving 
up into the domed area underneath the free surface if the calculations 
had proceeded further. Figure 13 shows a sample plot of the stream- 
lines, corresponding to the maximum time at which the calculations 
remained stable. It is important to note that the results cited above 
have not been subjected to any kind of numerical filtering. The cal- 
culations are generally terminated when the free surface forms a sharp 
corner or when difficulties are encountered in obtaining rapid conver- 
gence. There does not seem to be any direct relation between the 
breakdown of the numerical model and the attainment of maximum 
scar conditions. This is partly because the effect of viscous and turbu- 
lent diffusion on the counter-rotating vortices and on the scars and 
striations is not yet taken into account. Nevertheless, the idealized 
model gives some indication of the surface deformation and validates 
the conclusion that the rise of the vortex pair does not give rise to 
waves on either side of the scars, as anticipated by Kochin, et al. 
[Ref. 15]. Finally, it is clear that a two-dimensional numerical model 
cannot predict the existence of striations which are essentially three- 
dimensional instabilities. It is quite possible that the most important 
physical phenomenon leading to the existence and explanation of the 
dark narrow regions of wake behind the ships may be the striations or 


the interaction of striations with scars. 
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Ill. EXPERIMENTS AND COMPARISONS 


A. EXPERIMENTAL APPARATUS 


Experiments were conducted in a water basin which was twelve 


feet long, three feet wide, and four ieet deep. Additional details of the 


basin are described by Leeker [Ref. 22] and will not be repeated here. 


Two-dimensional trailing vortices were generated through the use 


of three different methods: 


|e 


a 


The use of two counter-rotating plates which, when fully rotated, 
were flush with the bottom of the tank, i.e., they have “vanished” 
at the end of their rotation, as far as the resulting fluid motion is 
concerned (see Figure 14), and no secondary vortices were 
generated. 


. The use of a two-dimensional, sharp-edged slot on the wall of a 


chamber, attached to the bottom of the basin. The fluid in the 
chamber is driven by a piston, attached to a pulley and weight 
system. 


The use of a streamlined nozzle attached to the aforementioned 
chamber, directly in front of the nozzle. Most of the experiments 
were conducted through the use of the nozzle because of the 
superior quality of the resulting vortices. 


To visualize the flow created by the vortices, a fluorescein dye was 


introduced either from the false bottom at about the midpoint of each 


rotating plate or from the inner walls of the nozzle. The dye was con- 


trolled by making it slightly buoyant with the use of alcohol and by 


using a flow control valve just outside the tank near the dye exit. 
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Figure 14. Rotating Plates and the Kelvin Oval 
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The data collection system consisted of a video camera, a VHS 
video cassette recorder, a monitor, and a date/timer. The date/timer 
was used to imprint on the screen and the tape a running clock and 


date (which was used as a run number) for each experiment. 


B. PROCEDURES 

The tank was first filled to a suitable level and allowed to sit until 
the fluid had come to rest. Water levels were varied throughout the 
experimentation to assess the effect of the water level on how soon 
the vortices became turbulent. While waiting for the fluid to come to 
rest, the plates were opened and the weights set. After the fluid was at 
rest, the dye was slowly introduced into the system. When the dye had 
risen just above the tips of the plates, the collimated light source was 
turned on and the video recording system was set to run. The weights 
were then actuated, causing the plates to rotate, which created the 
two-dimensional vortices. The vortices rose by their mutual induction 
velocity. As the vortices continued to rise, the free surface started to 


deform, leading to the formation of scars and striations. 


C. DATA ACQUISITION SYSTEM 

A NEC TI-22P CCD camera was used to record the experimental 
runs onto a Panasonic AG-6300 Video Cassette Recorder. The VHS 
recorder is required because for every frame of video information dis- 


played on the monitor, the VHS tape stores another frame of 
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information digitally which can later be extracted with the proper 
equipment. A Sun! computer with a motion analysis VP110 system was 
used to analyze the data. The motion analysis system allows the user to 
superimpose on the video monitor an outline of areas that have the 
same intensity or brightness levels. The brightness or area to be out- 
lined is controlled by the user as he sets a threshold. The threshold 
can be varied from zero (black) to 1999 (white). The area outlined 
shows up on the monitor, allowing the user fine control of what is to 
be picked up by the computer. 

The outlines are then turned into digital data and retained in a file 
by pixel location of the monitor screen. Once the data is stored, it can 
be edited and manipulated to eliminate unwanted images (see Figures 
15a and 15b). Further manipulation then allows the user to find cen- 
troids of areas, paths of centroids, and velocity and acceleration of the 
particles on these paths. This is all done through the use of the soft- 
ware provided in a package known as ExpertVision.? 

Figures 16 and 17 show the evolution of the free surface, as 
determined through the use of the ExpertVision, for the Froude num- 
bers of Fr = 0.39 and 0.53. Also noted in these figures are the posi- 
tions of the counter-rotating vortices. These will later be compared 


with those predicted numerically. 


1Sun is a registered trademark of Sun Microsystems, Inc. 


2ExpertVision is a registered trademark of Motion Analysis 
Corporation. 
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Figure 15a. Identification of the Free Surface and the Kelvin Oval 
(early stages) 
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Figure 15b. Identification of the Free Surface and the Kelvin Oval 
(later stages) 


36 





Figure 16. Evolution of the Free Surface During the Rise of 
Vortices for Fr = 0.39 
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Figure 17. Evolution of the Free Surface During the Rise of 
Vortices for Fr = 0.53 
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D. COMPARISON OF ANALYSIS AND EXPERIMENTS 

In general, the numerical model predicted closely several features 
that occurred in the experiments: the free surface shape, particularly 
the scar region and the smooth shape of the hump or domed area in 
the center. Furthermore, the model predicted that litthe would hap- 
pen beyond Y/bo = 3. This was found to be true for the stages of defor- 
mation up until the time when the vortices lost nearly all of their 
vertical movement. The model also predicted the location of the scars 
to a high degree of accuracy as well as the slight inward movement of 
the path of the vortices during their vertical ascent as they are drawn 
into the humped region. 

Figures 18 through 20 show a comparison of the variation of the 
local maximum of the free-surface deformation as a function of the 
instantaneous position of the counter rotating vortex pair for Fr = 
0.39, 0.42, and 0.53, respectively. Where applicable, the agreement 
between the experimental data and numerical predictions is quite 
satisfactory, particularly when one considers the fact that the surface 
deformations take place in a time period less than 0.3 seconds and the 
fact that the numerical model is an inviscid flow idealization of our 
conceptual view of the phenomenon. The smaller the number of vor- 
tices or marker elements, the larger the dome height and slower the 
onset of the numerical instability. This is neither surprising nor 
entirely unexpected. The stability of vortex sheets and, in particular, 


the Helmholtz instability have been the subject of intense interest. It is 
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Figure 18. Comparison of the Measured and Predicted Scar Elevations 
for Fr = 0.39 
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Figure 19. Comparison of the Measured and Predicted Scar 
Elevations for Fr = 0.42 
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Figure 20. Comparison of the Measured and Predicted Scar 
Elevations for Fr = 0.53 
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a well-known fact that an infinitesimal disturbance of wavelength A on 
a plane sheet of strength y grows like exp(zyt/A), according to which 
the shorter waves grow faster [Refs. 18, 23-25]. Clearly, the larger the 
number of vortices, the smaller the wave length and the faster and 
sooner the inception of instability, as demonstrated once again by the 
present calculations. As noted earlier, the instability could have been 
further delayed by filtering the numerical errors leading to the saw- 
tooth instability. In fact, this could have yielded scar elevations even 
larger than those observed experimentally. But the use of an artificial 
filter has been disregarded. The purpose of the present calculations 
was not the exact prediction of the ultimate scar shape but rather a 
search for new insights and inspiration in the simulation of a complex 
three-dimensional free-surface flow. The challenge was to predict the 
rate of growth in the face of ever-accumulating errors within a time 
scale appropriate for the investigation. The numerical model cannot 
and is not expected to deal with three-dimensional instabilities and 
turbulence seen during the actual experiments. Although the model 
has the capability to vary the vortex strength with time to account for 
the effects of laminar and turbulent diffusion, this was not done. By 
letting the vortex strength diminish with time, the program may have 
run longer and the actual height of the humped region may have 


agreed more closely with the heights seen during the experiments. 
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IV. CONCLUSIONS 


The investigation reported herein warranted the following 
conclusions: 


1. The numerical model, based on the use of a vortex sheet, redis- 
cretization of the sheet segments, and the desingularization of 
the complex velocity, without the use of an additional numerical 
filter, is capable of predicting the rate of growth of the free-sur- 
face shapes and the location of the early stages of the free surface 
deformations. The ultimate heights of the scars and striations are 
dictated by laminar and turbulent diffusion of vortices, surface 
impurities, and, to a lesser extent, by the initial shape of the 
counter-rotating vortex pair. The calculations are generally ter- 
minated when the free surface forms a sharp corner or when dif- 
ficulties are encountered in obtaining rapid convergence. There 
does not seem to be any direct relation between the conditions 
prevailing at the time of breakdown of the numerical model and 
the physical conditions existing at the time of attainment of 
maximum free-surface deformation. Clearly, there are other phe- 
nomena, such as turbulence and three-dimensional effects, with 
which the present numerical model cannot deal. 


2. For any of the Froude numbers encountered in the analysis or 
experiments, no wave train was observed on either side of the 
scars. This conclusion was somewhat anticipated on the basis of 
previous calculations. For very small Froude numbers (say, 
Fr < 0.15), the scars were very small, the vortices followed the 
rigid-wall path, and the scars were slaved to the vortices. 


3. Preliminary experiments have revealed the existence of striations 
even in an otherwise two-dimensional vortex flow. In other 
words, the striations observed in earlier investigations at the 
Naval Postgraduate School with an inclined trailing vortex pair 
[Ref. 26] were not due to the three-dimensionality of the flow. It 
appears that the striations are a consequence of gravitational and 
centrifugal instability, in both two- and three-dimensional vortical 
flows occurring in the vicinity of a free surface. It also appears 
that the turbulent wake-like features observed in SAR pictures 
are related to a vortex mechanism: formation of scars and Stria- 
tions, intersection of the striation waves with the edges of the 
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scars, and the formation of small-scale vortical patterns. These 
features will be the subject of continuing investigation. 
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